JOURNAL  OF  GEOPHYSICAL  RESEARCH,  VOL.  ???,  XXXX,  DOL10.1029/, 


Decadal  Variability  of  the  Tropical  Stratosphere: 
Secondary  Influence  of  the  El  Nino-Southern 
Oscillation 

L.  L.  Hood  and  B.  E.  Soukharev 

Lunar  and  Planetary  Laboratory,  University  of  Arizona,  Tucson,  Arizona, 

USA 

J.  P.  McCormack 

Naval  Research  Laboratory,  Washington,  D.  C.,  USA 


L.  L.  Hood  and  B.  E.  Soukharev,  Lunar  and  Planetary  Laboratory,  University 
of  Arizona,  1629  E.  University  Blvd.,  Tucson,  AZ  85721,  USA  (lon@lpl.arizona.edu; 
boris_soukharev@hotmail.com);  J.  P.  McCormack,  Naval  Research  Laboratory,  Code  7646,  4555 
Overlook  Ave.  SW,  Washington,  D.  C.,  20375,  USA  (mccormack@nrl.navy.mil). 


DRAFT 


February  4,  2010,  10:02am 


DRAFT 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

2010 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-2010  to  00-00-2010 


4.  TITLE  AND  SUBTITLE 

Decadal  Variability  of  the  Tropical  Stratosphere:  Secondary  Influence  of 
the  El  Nino-Southern  Oscillation 

6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory, Code  7646,4555  Overlook  Ave. 
SW, Washington, DC, 20375 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROIECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

Journal  Geophysical  Research,  in  press,  2010 

14.  ABSTRACT 

see  report 

15.  SUBIECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OF 

18.  NUMBER 

ABSTRACT 

OF  PAGES 

Same  as 

61 

Report  (SAR) 

19a.  NAME  OF 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


X  -  2 


HOOD  ET  AL.:  DECADAL  VARIABILITY  OF  THE  TROPICAL  STRATOSPHERE 


Abstract.  A  decadal  variation  of  tropical  lower  stratospheric  ozone  and 
temperature  has  previously  been  identified  that  correlates  positively  with  the 
11-year  solar  activity  cycle.  However,  the  El  Nino-Southern  Oscillation  (ENSO) 
also  influences  lower  stratospheric  ozone  /  temperature.  It  is  therefore  legit¬ 
imate  to  ask  whether  quasi-decadal  ENSO  variability  can  contribute  to  this 
apparent  solar  cycle  variation,  either  accidentally  because  of  the  short  mea¬ 
surement  record  (e.g.,  Marsh  and  Garcia  [2007])  or  physically  because  so¬ 
lar  variability  affects  ENSO.  Here,  we  present  multiple  regression  analyses 
of  available  data  records  to  compare  differences  in  results  obtained  with  and 
without  including  an  ENSO  term  in  the  statistical  model.  In  addition,  sim¬ 
ulations  are  performed  using  the  NRL  NOGAPS-ALPHA  GCM  for  warm/cold 
ENSO  conditions  to  test  for  consistency  with  the  ENSO  regression  results. 

We  find  only  very  minor  changes  in  annual  mean  solar  regression  coefficients 
when  an  ENSO  term  is  included.  However,  the  observed  tropical  ENSO  re¬ 
sponse  provides  useful  insights  into  the  origin  of  the  unexpected  vertical  struc¬ 
ture  of  the  tropical  solar  cycle  ozone  response.  In  particular,  the  ENSO  ozone 
response  is  negative  in  the  lower  stratosphere  due  to  increased  upwelling  but 
changes  sign,  becoming  positive  in  the  middle  stratosphere  (5  —  10  hPa) 
due  mainly  to  advective  decreases  of  temperature  and  NOx,  which  photo- 
chemically  increase  ozone.  A  similar  mechanism  may  explain  the  observed 
lower  stratospheric  solar  cycle  ozone  /  temperature  response  and  the  absence 
of  a  significant  response  in  the  tropical  middle  stratosphere. 
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1.  Introduction 

Continuous  global  satellite  measurements  beginning  in  late  1978  show  a  decadal  vari¬ 
ation  of  total  column  ozone  at  tropical  and  subtropical  latitudes  that  is  approximately 
in  phase  with  the  solar  cycle  (Figure  lb).  Observations  also  indicate  a  significant  solar 
cycle  temperature  variation  in  both  the  upper  stratosphere  [e.g.,  McCormack  and  Hood, 
1996]  and  in  the  lower  stratosphere  [e.g.,  Steinbrecht  et  ah,  2003;  Labitzke,  2004;  Crooks 
and  Gray,  2005;  Randel  et  ah,  2009a;  Gray  et  ah,  2009].  However,  the  shortness  of  the 
record  and  the  occurrence  of  two  major  volcanic  eruptions  following  the  first  two  solar 
maxima  (El  Chichon  in  early  April  of  1982  and  Pinatubo  in  June  of  1991)  has  led  to 
questions  about  whether  the  ozone  response  to  volcanic  aerosol  injections  could  alias  or 
even  be  confused  with  the  response  to  the  11-year  solar  cycle  [Solomon  et  ah,  1996;  Lee 
and  Smith,  2003] .  This  is  especially  true  since  the  statistical  method  that  has  been  most 
commonly  applied  to  estimate  empirically  the  solar  component  of  stratospheric  interan¬ 
nual  variability  is  multiple  linear  regression,  which  makes  the  simplified  assumption  that 
all  forcings  are  orthogonal,  i.e.,  independent  of  one  another,  and  that  all  responses  are 
linear  [e.g.,  Stolarski  et  ah,  1991;  Hood  and  McCormack,  1992;  Hood  et  ah,  1993;  Reinscl 
et  ah,  2005]. 

These  concerns  have  fortunately  been  mitigated  in  recent  years  by  the  absence  of  major 
eruptions  since  1991  and  a  continuation  of  the  decadal  ozone  variation  for  one  additional 
cycle.  Moreover,  extensions  of  the  total  ozone  record  backward  in  time  using  ground-based 
Dobson  spectrophotometer  data  have  supported  the  existence  of  a  decadal  ozone  variation 
for  at  least  several  more  cycles  prior  to  the  start  of  continuous  satellite  data  [WMO,  2007, 
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Chapter  3].  Consequently,  the  existence  of  a  significant  solar  cycle  variation  of  total  ozone 
(and,  by  extension,  the  stratosphere  as  a  whole)  has  become  clearer  in  recent  years  [e.g., 
WMO,  2007;  Gray  et  al.,  2010], 

Nevertheless,  the  vertical  structure  of  the  stratospheric  solar  cycle  variation  remains 
poorly  understood.  Prior  to  2006,  most  or  all  published  model  simulations  of  the  tropical 
solar  cycle  ozone  variation  were  characterized  by  a  broad  response  maximum  centered  in 
the  middle  stratosphere  (5  -  10  hPa)  that  was  caused  mainly  by  increased  photolysis  of 
molecular  oxygen  by  solar  UV  radiation  (see,  e.g.,  Hood  [2004]  and  Chapter  3  of  WMO 
[2007]  for  reviews).  In  contrast,  the  observationally  estimated  tropical  ozone  response  is 
characterized  by  a  double-peaked  vertical  structure  with  solar  min-to-max  amplitudes  of 
several  per  cent  in  the  upper  and  lower  stratosphere  but  no  statistically  significant  response 
in  the  middle  stratosphere  (e.g.,  McCormack  and  Hood  [1996]).  It  was  not  clear  whether 
the  disagreement  between  models  and  observations  was  due  to  a  model  shortcoming  or  to 
problems  with  the  observations  [e.g.,  Lee  and  Smith,  2003]. 

To  address  this  issue,  an  improved  analysis  of  three  independent  satellite  ozone  datasets 
with  lengths  extending  up  to  25  years  was  carried  out  [Soukharev  and  Hood,  2006  (here¬ 
after  SH06);  see  also  Randel  and  Wu,  2007].  Column  ozone  measurements  were  also 
compared  with  LIpper  Atmosphere  Research  Satellite  (UARS)  Halogen  Occultation  Ex¬ 
periment  (HALOE)  ozone  profile  data  during  the  1992-2003  period  when  no  major  vol¬ 
canic  eruptions  occurred  to  test  the  conclusion  that  most  of  the  solar  cycle  total  ozone 
variation  occurred  in  the  lower  stratosphere.  Separate  time  intervals  were  analyzed  us¬ 
ing  the  longest  (25-year)  data  record  to  test  the  reproducibility  of  the  solar  regression 
coefficients.  Results  conhrmed  that  the  basic  vertical  structure  of  the  derived  solar  cycle 
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ozone  response  was  the  same  during  the  last  half  (post-Pinatubo)  part  of  the  record  as  it 
was  during  the  first  half  and  that  most  of  the  solar  cycle  column  ozone  variation  occurred 
in  the  lower  stratosphere.  Overall,  therefore,  the  analysis  supported  the  reality  of  the 
double-peaked  ozone  profile  response  to  the  solar  cycle. 

Consistent  with  the  weight  of  the  observational  evidence,  simulations  by  coupled  chem¬ 
istry  climate  models  (CCM’s)  have  begun  to  produce  a  double-peaked  ozone  profile  re¬ 
sponse  that  is  qualitatively  similar  to  that  which  is  observed  [e.g.,  Austin  et  ah,  2008; 
Matthes  et  ah,  2007].  The  main  improvements  in  the  latter  simulations  were:  (a)  In¬ 
tegration  of  each  CCM  forward  in  time  (“transient  simulations”);  and  (b)  forcing  the 
models  with  observed  sea  surface  temperatures  (SST’s).  However,  the  specific  physical 
mechanisms  that  are  responsible  for  this  better  agreement  are  difficult  to  identify  from 
the  CCM  simulations  alone. 

Most  current  mechanisms  for  explaining  a  solar  cycle  variation  in  the  lower  stratosphere 
involve  “downward  control”  from  the  upper  atmosphere.  These  include  (a)  influences  of 
observed  11-year  variations  in  the  lower  mesospheric  subtropical  jet  on  the  timing  (early 
winter  vs.  late  winter)  of  sudden  stratospheric  warmings,  which  in  turn  modulate  the 
Brewer-Dobson  circulation  [Kodera  and  Kuroda,  2002;  Matthes  et  al.,  2004;  2006;  Gray  et 
al.,  2006;  Ito  et  ah,  2009];  and  (b)  solar-induced  changes  in  the  duration  of  the  west  phase 
of  the  QBO  [McCormack  et  ah,  2007,  and  references  therein].  However,  it  is  also  possible 
that  “feedbacks  from  below”  are  a  significant  contributing  cause.  Recent  experiments 
using  global  climate  models  have  shown  that  both  solar  UV  forcing  of  the  stratosphere  and 
an  amplifying,  coupled  ocean-troposphere  response  (driven  by  indirect  effects  of  solar  UV 
forcing  as  well  as  direct  effects  of  total  irradiance  forcing)  must  be  included  to  simulate 
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the  magnitude  of  observed  11-year  climate  signals  in  the  tropical  Pacific  [Meehl  et  ah, 
2009].  The  observationally  estimated  surface  response  to  positive  maxima  in  11-year  solar 
forcing  resembles  that  observed  during  cold  ENSO  events  although  there  are  important 
differences  [van  Loon  et  ah,  2007;  van  Loon  and  Meehl,  2008].  Other  authors  have  also 
obtained  observational  evidence  for  a  small  solar  cycle  signal  in  global  SST’s  [White  and 
Liu,  2008,  and  references  therein]  as  well  as  in  equatorial  tropospheric  temperatures  [Salby 
and  Callaghan,  2006].  In  contrast,  there  is  no  evidence  for  solar  forcing  of  ENSO  itself 
since  solar  flux  and  ENSO  indices  do  not  correlate  significantly  over  long  (>  25-year) 
time  periods.  Nevertheless,  these  observational  and  climate  model  studies  suggest  that 
feedbacks  from  the  amplifying  ocean-troposphere  response  may  be  involved  in  producing 
the  unexpected  lower  stratospheric  ozone  and  temperature  responses.  This  possibility 
may  be  consistent  with  the  need  to  use  observed  SST’s  as  a  lower  boundary  condition 
in  CCM  simulations  that  successfully  produce  a  second  lower  stratospheric  solar  ozone 
response  maximum  [e.g.,  Austin  et  al.,  2008]. 

However,  a  study  using  the  National  Center  for  Atmospheric  Research  (NCAR)  Whole 
Atmosphere  Community  Climate  Model  (WACCM)  has  suggested  that  aliasing  from  the 
El-Nino  /  Southern  Oscillation  (ENSO)  could  be  causing  part  of  the  apparent  positive 
solar  cycle  ozone  response  in  the  tropical  lower  stratosphere  [Marsh  and  Garcia,  2007].  It 
was  pointed  out  that,  over  the  1979  -  2003  time  period,  the  Nino  3.4  ENSO  index  (N3.4) 
was  almost  significantly  correlated  (R  ~  -  0.19  at  a  lag  of  —6  months)  with  solar  UV 
proxies  such  as  the  10.7  cm  radio  flux  (F10.7).  Since  F10.7  and  N3.4  may  not  therefore 
have  been  entirely  orthogonal  during  this  time  period,  some  aliasing  of  the  ENSO  response 
onto  the  solar  cycle  response  could  have  occurred.  In  support  of  this  possibility,  multiple 
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regression  analyses  of  WACOM  model  output  with  and  without  including  an  ENSO  proxy 
term  with  a  suitable  time  lag  in  the  regression  model  showed  that  including  the  ENSO 
term  reduced  considerably  the  apparent  solar  cycle  ozone  response  at  levels  below  20  hPa. 
Since  the  statistical  model  of  SH06  did  not  include  an  ENSO  term,  the  same  aliasing 
could  potentially  exist  for  the  ozone  response  derived  from  observations. 

In  this  paper,  new  multiple  regression  analyses  are  reported  of  available  long-term  satel¬ 
lite  remote  sensing  data  records  using  an  improved  statistical  model  that  includes  an 
ENSO  term.  In  addition,  time-progressive  simulations  using  the  Naval  Research  Labora¬ 
tory  (NRL)  NOGAPS-ALPHA  general  circulation  model  (GCM)  are  carried  out  for  warm 
and  cold  ENSO  conditions  to  allow  comparisons  with  the  ENSO  regression  results.  The 
main  objectives  are  (a)  to  investigate  whether  aliasing  from  ENSO  contributes  signifi¬ 
cantly  to  the  solar  cycle  ozone  response  derived  from  observations;  and  (b)  to  investigate 
more  completely  the  response  of  the  tropical  stratosphere  to  ENSO  forcing  using  both  ob¬ 
servational  analyses  and  model  simulations.  We  find  that  the  ENSO  response  in  both  the 
lower  and  middle  stratosphere  provides  useful  insights  into  the  origin  of  the  unexpected 
vertical  structure  of  the  tropical  solar  cycle  ozone  response. 

The  paper  is  organized  as  follows.  In  section  2,  the  existence  of  both  solar  cycle  and 
ENSO  components  of  decadal  stratospheric  variability  is  first  illustrated  using  representa¬ 
tive  total  ozone  time  series.  The  adopted  multiple  regression  statistical  model  is  described 
and  the  existence  of  any  significant  aliasing  of  the  solar  cycle  ozone  profile  response  by 
ENSO  is  tested  by  analyzing  data  with  and  without  inclusion  of  an  ENSO  term.  In  sec¬ 
tion  3,  the  statistical  model  is  applied  to  estimate  ENSO  regression  coefficients  for  the 
ozone  profile,  temperature,  zonal  wind,  and  total  ozone  using  available  long-term  data 
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sets.  The  geographic  and  seasonal  dependences  of  the  total  ozone  ENSO  coefficient  are 
also  estimated.  In  section  4,  model  calculations  for  strong  and  weak  ENSO  conditions  are 
carried  out  using  the  NOGAPS-ALPHA  model  to  allow  comparisons  with  the  observa- 
tionally  derived  ENSO  regression  coefficients.  In  section  5,  likely  physical  mechanisms  for 
explaining  the  vertical  structures  of  the  ENSO  regression  results  are  discussed.  Possible 
implications  for  the  origin  of  the  vertical  structure  of  the  solar  cycle  ozone  response  are 
also  examined.  A  summary  and  further  discussion  are  given  in  section  6. 


2.  Statistical  Separation  of  the  ENSO  and  Solar  Responses 

Figure  1  compares  several  zonal  averages  of  the  Total  Ozone  Mapping  Spec¬ 
trometer  /  Solar  Backscattered  Ultraviolet  (TOMS/SBUV)  total  ozone  record  com¬ 
piled  and  calibrated  at  Goddard  Space  Flight  Center  by  S.  Frith  and  R.  Stolarski 
(http://code916.gsfc.nasa.gov/Data_services;  Frith  and  Stolarski  [2005])  with  the  solar 
Mg  II  UV  index  (ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SOLAR_UV;  Viereck  and 
Puga  [1999]).  In  Figure  la,  which  plots  a  nearly  global  (65°S  to  65°N)  average  of  the 
data,  at  least  three  sources  of  decadal  variability  can  be  seen.  First,  there  is  a  long-term, 
at  least  partly  anthropogenic,  trend.  This  trend  was  roughly  linear  and  negative  until 
the  late  1990’s  but  has  since  been  noticeably  reduced  in  amplitude  [WMO,  2007,  Chapter 
3].  Second,  significant  ozone  decreases  occur  following  the  major  volcanic  eruptions  of  El 
Chichon  in  late  March  /  early  April  of  1982  and  Pinatubo  in  June  of  1991  [Solomon  et  ah, 
1996;  Robock,  2000;  Stenchikov  et  al.,  2002],  Third,  there  is  a  tendency  for  higher  global 
ozone  averages  near  maxima  of  the  11-year  solar  cycle  [Hood  and  McCormack,  1992;  Hood, 
1997;  Zerefos  et  ah,  1997;  Soukharev  and  Hood,  2006;  Randel  and  Wu,  2007;  Tourpali  et 
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al.,  2007].  The  volcanic  and  trend  contributions  are  strongest  at  middle  and  high  latitudes 
where  heterogeneous  chemical  losses  and  dynamical  transport  are  most  important  (e.g., 
Brasseur  and  Solomon  [2005]).  However,  at  low  latitudes,  the  solar  contribution  appears 
to  dominate.  This  is  illustrated  in  Figure  lb,  which  plots  a  tropical  and  subtropical  (35°S 
to  35°N)  average  time  series. 

The  El-Nino  /  Southern  Oscillation  is  an  oscillation  of  the  coupled  atmosphere  and  ocean 
system  (primarily  in  the  tropical  Pacific)  that  is  arguably  the  dominant  global  mode  of 
interannual  climate  variability  [e.g.,  Bronnimann,  2007].  Warmer  SST’s  during  an  El-Nino 
event  produce  increased  convection,  increased  latent  heat  release,  a  warmer  troposphere, 
and  a  higher  tropopause  particularly  over  the  eastern  tropical  Pacific  but  also  throughout 
the  tropics.  The  opposite  is  the  case  for  La  Nina.  Evidence  for  an  ENSO  contribution  to 
total  ozone  variability  in  the  tropics  was  first  obtained  using  Nimbus  4  BUV  data  [Hasebe, 
1983]  and  the  geographic  dependence  of  the  variation  was  first  estimated  using  Nimbus  7 
TOMS  data  [Shiotani,  1992;  Randel  and  Cobb,  1994], 

Figure  2  compares  a  tropical  zonal  average  (20°S  to  20°N)  of  the  Version  8  TOMS/SBLTV 
record  with  the  Nino  3.4  index  (sea  surface  temperature  anomaly  averaged  over  the  5°S  to 
5°N  and  120°W  to  170°W  sector).  A  number  of  distinct  positive  (El-Nino)  and  negative 
(La  Nina)  events  can  be  seen  in  Figure  2b.  Note  that  El-Nino  peaks  tend  to  occur  in 
early  boreal  winter.  The  peaks  in  the  N3.4  index  in  Figure  2b  are  often  accompanied 
by  decreases  several  months  later  in  ozone  in  Figure  2a.  For  example,  a  sharp  tropical 
ozone  decrease  occurs  several  months  following  the  El-Nino  event  of  late  1997.  This 
ozone  reduction  is  generally  understood  to  be  caused  by  increased  tropical  upwelling 
following  warm  events  (see  section  5.1  below).  However,  dynamically  induced  interannual 
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ozone  variability  from  other  sources  (notably  from  the  equatorial  quasi-biennial  wind 
oscillation,  or  QBO)  is  large  in  the  tropics.  Consequently,  the  best  correlation  between 
tropical  TOMS/SBUV  zonal  mean  column  ozone  and  the  Nino  3.4  ENSO  index  is  only 
R  =  —0.37  (significant  at  99%  confidence)  at  4  months  positive  lag  (ozone  lags  ENSO). 
In  general,  Figures  1  and  2  suggest  that  ENSO  is  secondary  in  importance  to  solar  and 
volcanic  forcing  for  driving  decadal  variability  of  the  tropical  stratosphere. 

2.1.  Regression  Model 

In  order  to  represent  the  temporal  behavior  of  a  given  stratospheric  variable,  X(t),  we 
consider  a  multiple  linear  regression  model  with  a  form  similar  to  that  used  in  previous 
work  (SH06  and  references  therein),  but  with  the  addition  of  an  ENSO  term  and  including 
a  more  complete  representation  of  the  QBO: 

%  (t)  =  H\i)  +  fltrendt  +  (^QBOlu30hPa{t  —  ^QROl)  +  f^QB02UlQhPa{t  —  lagQB02) 

+ PvoicanicAer  osol{t)  +  (5soiarM  gl  I  (t)  +  (3ENSoN3A{t  —  lagENSO )  +  e{t)  (1) 

where  t  is  the  time  in  months  (for  ERA-40  and  TOMS/SBUV)  or  3- month  seasonal 
increments  (for  ozone  profile  datasets),  /i(i)  is  the  long-term  mean  for  the  ith  month 
( i  =  1,2,. ..12)  or  season  ( i  =  1,2,  ...4)  of  the  year;  u^hPa  and  uiohPa  are  the  30 
and  10  hPa  equatorial  zonal  winds  obtained  from  National  Centers  for  Environmen¬ 
tal  Prediction  (NCEP)  reanalysis  data  [Kalnay  et  al.,  1996];  lagQBoi  and  lagQB02  are 
lag  times  required  to  produce  a  maximum  positive  or  negative  correlation  between  the 
X 1(f)  time  series  at  a  given  location  and  the  30  and  10  hPa  equatorial  winds;  Aerosol{t ) 
is  a  stratospheric  aerosol  index  based  on  a  combination  of  SAM  II  and  Stratospheric 
Aerosol  and  Gas  Experiment  (SAGE)  I /II  optical  depth  measurements  [Thomason  et 
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al.,  1997];  MgII(t)  is  the  core-to-wing  ratio  of  the  solar  Mg  II  line  at  280  nm  based 
on  Nimbus  7  SBUV,  NOAA  9  and  11  SBUV/2,  and  UARS  SUSIM  data  [e.g.,  Viereck 
and  Puga,  1999];  N3.4  is  the  Nino  3.4  index  (defined  as  above  and  available  from: 
http://www.cgd.ucar.edu/cas/catalog/climind/Nino_3_3.4_indices.html);  lagENSO  is  a  lag 
time  required  to  produce  a  maximum  correlation  amplitude  between  X(t)  and  N3.4;  and 
e(t)  is  a  residual  error  term.  The  coefficients  /3trend,  Pqboi ,  Pqbo2 ,  Pvoicanic ,  Psoiar,  and 
Pen  so  are  determined  by  least  squares  regression.  The  aerosol  index  term  is  included  in 
the  model  only  for  levels  at  10  hPa  and  below. 

As  in  previous  work,  to  avoid  overestimation  of  regression  coefficients  and  their  statis¬ 
tical  significance,  it  is  important  to  model  the  residual  error  term  as  an  autoregressive 
process,  e.g.,  for  a  first-order  process,  e{t)  =  re(t  —  1)  +  w(t),  where  w(t)  is  white  noise 
and  r  is  determined  in  an  initial  application  of  (1)  (e.g.,  Neter  et  al.  [1985]).  Crooks 
and  Gray  [2005]  have  modeled  the  error  term  as  a  third-order  autoregressive  process  to 
be  more  complete;  however,  they  find  no  difference  in  the  results  if  other  orders  are  used 
instead.  A  first-order  process  is  therefore  considered  to  be  sufficient  here. 

There  are  several  alternate  methods  for  representing  the  QBO  in  multiple  regression 
statistical  models.  A  single  QBO  term  proportional  to  the  30  hPa  equatorial  wind  and 
with  an  optimum  phase  lag  has  been  found  to  be  adequate  for  evaluating  the  solar  cycle 
ozone  variation  [SH06].  However,  for  the  purpose  of  evaluating  solar  cycle  or  ENSO  zonal 
wind  and  temperature  variations,  a  more  complete  QBO  representation  is  necessary.  In 
particular,  it  is  preferable  to  include  two  separate  QBO  terms  that  account  for  the  out-of- 
phase  relationship  between  QBO  winds  in  the  upper  and  lower  stratosphere  (e.g.,  Randel 
and  Wu  [1996];  Crooks  and  Gray  [2005];  Shibata  and  Deushi  [2008]).  Randel  and  Wu 
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[2007;  1996]  have  used  the  first  two  QBO  empirical  orthogonal  functions  (EOF’s)  derived 
from  equatorial  zonal  wind  data  at  all  levels  to  represent  the  QBO  in  their  regression 
model.  On  the  other  hand,  Shibata  and  Deushi  [2008]  and  Steinbrecht  et  al.  [2003] 
have  used  the  actual  QBO  equatorial  winds  at  two  different  levels  to  represent  the  QBO. 
The  former  authors  used  50  and  20  hPa  winds  while  the  latter  authors  used  30  and  10 
hPa  winds.  The  winds  at  either  of  these  two  levels  are  about  one  quarter  cycle  out  of 
phase  indicating  that  they  can  be  treated  as  approximately  independent  variables.  In  the 
present  work,  we  follow  Steinbrecht  et  al.  and  use  actual  QBO  equatorial  winds  at  the 
30  and  10  hPa  levels  but  with  small  lag  times  added  to  further  maximize  the  projection 
of  X(t )  onto  the  QBO.  Several  experimental  analyses  were  carried  out  using  either  the 
orthogonal  time  series  of  Randel  and  Wu  (private  communication,  2008)  or  the  actual 
winds  at  various  levels  (50  and  20  hPa,  30  and  10  hPa).  It  was  found  that  there  is  no 
statistically  significant  difference  in  the  results  for  the  solar  and  ENSO  zonal  wind  and 
temperature  regression  coefficients  when  any  of  these  three  methods  is  applied. 

2.2.  Tests  for  ENSO  Aliasing 

The  zonal  mean  ozone  profile  response  to  the  11-year  solar  activity  cycle  has  previously 
been  estimated  by  a  number  of  analysts  using  mainly  satellite  remote  sensing  data  (e.g., 
Hood  et  al.  [1993];  Chandra  and  McPeters  [1994];  McCormack  and  Hood  [1996];  Wang 
et  al.  [1996];  Soukharev  and  Hood  [2006];  Randel  and  Wu  [2007];  Tourpali  et  al.  [2007]). 
However,  as  mentioned  in  the  Introduction,  it  has  been  found  using  WACCM  model 
output  that  the  solar  cycle  ozone  regression  coefficient  in  the  lower  stratosphere  decreases 
significantly  when  an  ENSO  term  is  included  in  the  statistical  model,  suggesting  aliasing 
between  the  ENSO  and  solar  cycle  terms.  To  test  whether  this  also  occurs  in  actual 
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stratospheric  data,  we  have  re-calculated  the  annual  mean  solar  regression  coefficient 
using  the  same  data  sets  analyzed  previously  by  SH06. 

Figure  3  shows  results  of  an  application  of  (1)  with  and  without  inclusion  of  the  ENSO 
term  to  the  3- month  SBUV(/2)  zonal  mean  ozone  profile  time  series  over  the  1979  to 
2003  period  (see  SH06  for  a  description  of  this  data  set).  The  annual  mean  regression 
coefficient  obtained  without  the  ENSO  term  (top  panel)  is  compared  to  that  obtained 
with  the  ENSO  term  at  an  optimum  lag  of  +3  months  (bottom  panel).  They  are  nearly 
identical  and  there  is  no  significant  change  in  the  coefficients  in  the  lower  stratosphere. 
The  vertical  structure  in  the  tropics,  consisting  of  positive  responses  in  the  upper  and 
lower  stratosphere  but  no  signihcant  response  in  the  middle  stratosphere,  differs  from 
that  expected  from  photochemical  effects  alone  [SH06];  possible  explanations  for  this 
difference  are  discussed  in  section  5.3.  Similar  applications  to  other  stratospheric  data 
(SAGE  11  and  HALOE  ozone  profile  data,  ECMWF  ERA-40  temperatures  and  winds) 
also  yield  only  slightly  altered  annual  mean  solar  regression  coefficients  with  and  without 
inclusion  of  an  ENSO  term  (not  shown  here).  Figure  4  directly  compares  the  tropically 
averaged  solar  cycle  ozone  regression  coefficients  calculated  from  the  three  ozone  satellite 
data  sets  (SBUV,  UARS  HALOE,  and  SAGE  II)  using  (1)  with  and  without  an  ENSO 
term  at  optimum  lag  (3  months  since  the  time  series  are  seasonal  rather  than  monthly). 
The  changes  that  result  from  addition  of  the  ENSO  term  are  all  much  less  than  the  2 a 
error  bars.  Comparisons  of  the  SBUV(/2)  results  with  those  of  SAGE  II  and  HALOE 
are  necessary  since  SBLIV  does  not  adequately  resolve  the  ozone  profile  in  the  lowermost 
stratosphere  [Bhartia  et  ah,  1996;  SH06]. 
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We  have  also  investigated  whether  any  significant  aliasing  between  the  ENSO  and  solar 
cycle  terms  could  occur  when  only  data  for  individual  seasons  are  analyzed.  For  three 
of  the  four  seasons,  solar  cycle  ozone  regression  coefficients  obtained  with  and  without 
an  ENSO  term  are  nearly  identical,  as  is  the  case  in  Figures  3  and  4.  However,  for 
northern  winter  data,  some  significant  differences  are  found  (Figure  5).  Specifically,  the 
DJF  responses  in  the  lower  stratosphere  where  (1)  is  applied  with  an  ENSO  term  (Figure 
5b)  are  weakened  between  ~  20°S  and  30°N  by  as  much  as  25%  and  the  regions  of 
statistical  significance  are  smaller.  Also,  in  the  upper  stratosphere,  the  maximum  solar 
cycle  ozone  response  increases  from  ~  4%  to  ~  5%.  As  mentioned  above  in  relation  to 
Figure  2,  ENSO  events  tend  to  peak  during  boreal  winter  and  have  especially  large  effects 
in  the  eastern  Pacific  region.  It  is  therefore  not  unexpected  that  aliasing  from  ENSO 
would  preferentially  occur  during  this  season  and/or  at  certain  longitudes  where  ENSO 
forcing  is  more  prominent. 

As  a  further  test,  column  ozone  solar  regression  coefficients  were  calculated  as  a  function 
of  longitude,  latitude,  and  season  using  the  version  8  TOMS/SBUV  data  set.  Again,  for 
the  boreal  spring,  summer,  and  fall  seasons,  only  very  minor  differences  were  noted  when 
an  ENSO  term  was  added  to  the  regression  model.  However,  during  boreal  winter  (DJF), 
more  significant  differences  were  noted.  As  seen  in  Figure  6,  the  resulting  regression 
coefficients  are  slightly  altered  and  the  total  region  where  the  coefficients  are  significant 
(shaded  area)  is  somewhat  larger  when  the  ENSO  term  is  included.  However,  the  general 
character  and  distribution  of  the  DJF  solar  regression  coefficient  are  approximately  the 
same  in  Figures  6a  and  6b. 
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On  this  basis,  it  is  concluded  that  aliasing  from  ENSO  produces  only  minor  changes  in 
annual  mean  solar  cycle  regression  coefficients  calculated  from  available  stratospheric  data 
during  the  1979-2003  period.  The  annual  mean  and  most  of  the  seasonal  ozone  profile  solar 
regression  coefficients  obtained  previously  (e.g.,  SH06)  are  therefore  still  approximately 
valid.  Only  during  boreal  winter  (DJF)  are  significant  evidences  obtained  for  limited  (< 
25%)  aliasing  of  solar  regression  coefficients  at  certain  latitudes,  longitudes,  and  pressure 
levels.  The  overall  lack  of  sensitivity  of  the  observational  solar  cycle  regression  coefficients 
to  whether  an  ENSO  term  is  included  in  the  statistical  model  is  partly  because  of  other 
sources  of  interannual  variability  (e.g.,  the  QBO)  in  actual  stratospheric  time  series,  which 
obscure  the  ENSO  influence.  Another  contributing  factor  is  probably  that  the  WACOM 
model  overestimates  somewhat  the  actual  ozone  sensitivity  to  ENSO  in  parts  of  the  lower 
stratosphere  (see  next  section). 


3.  ENSO  Regression  Coefficients 

Figure  7  plots  the  annual  mean  Version  8  SBUV(/2)  ozone  profile  ENSO  regression 
coefficient  for  the  1979-2003  period,  expressed  in  per  cent  change  of  ozone  for  a  1  unit 
(i.e.,  1°  C)  positive  change  of  the  N3.4  index.  The  plot  was  produced  by  applying  (1) 
successively  to  the  3-month  average  SBUV(/2)  data  at  individual  latitude/pressure  grid 
points.  It  was  found  empirically  that  a  maximum  correlation  between  the  N3.4  index 
and  the  ozone  time  series  is  obtained  at  a  majority  of  grid  points  for  a  positive  phase  lag 
(Lenso)  of  1  time  unit  (3  months).  The  responses  shown  are  therefore  at  a  phase  lag  of  3 
months  (ozone  lags  ENSO).  Shaded  areas  indicate  grid  points  where  the  ENSO  regression 
coefficient  is  significant  at  the  95%  confidence  level.  As  seen  in  Figure  2,  a  1  unit  increase 
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in  N3.4  corresponds  to  a  modest  warm  (El-Nino)  event  in  the  eastern  tropical  Pacific  while 
a  2  to  2.5  unit  increase  corresponds  to  a  strong  El-Nino  event.  Therefore,  to  estimate  the 
true  ozone  response  to  a  strong  El-Nino  event,  the  values  in  the  tropical  lower  and  middle 
stratosphere  shown  in  Figure  7  should  be  multiplied  by  a  factor  of  2  or  2.5.  Figure  8  plots 
the  corresponding  regression  coefficients  obtained  using  SAGE  II  data  over  the  1985-2003 
period  and  UARS  HALOE  data  over  the  1992-2003  period.  For  the  SAGE  II  analysis, 
the  period  immediately  following  the  Pinatubo  eruption  (from  JJA  1991  to  SON  1993)  is 
excluded  from  consideration  in  lieu  of  including  an  aerosol  term  in  the  regression  model 
[SH06] . 

In  Figure  7,  a  large  region  of  significant  negative  SBUV(/2)  response  is  obtained  in  the 
tropical  lower  stratosphere  while  a  zone  of  positive  response  is  obtained  in  the  equatorial 
middle  stratosphere.  Similar  patterns  are  seen  in  Figures  8a  and  8b,  although  less  distinct, 
apparently  due  to  the  reduced  record  lengths  and/or  sampling  of  the  SAGE  II  and  HALOE 
data.  In  all  three  plots,  but  especially  in  the  SAGE  II  and  HALOE  plots,  the  zone  of 
positive  response  tends  to  be  shifted  noticeably  toward  the  Northern  Hemisphere.  As 
will  be  discussed  in  section  5,  both  the  negative  ozone  response  in  the  lower  stratosphere 
and  the  positive  response  in  the  equatorial  middle  stratosphere  are  consistent  with  a  net 
positive  increase  in  tropical  upwelling  rate  following  El-Nino  events.  In  addition,  a  zone  of 
apparently  significant  positive  response  is  obtained  at  low  latitudes  above  the  stratopause 
in  Figure  7.  However,  no  corresponding  response  is  seen  in  the  SAGE  II  and  HALOE 
plots  of  Figure  8. 

Table  1  lists  the  tropical  mean  (24°S  to  24°N)  ENSO  ozone  regression  coefficients  at 
five  lower  stratospheric  levels  as  calculated  here  by  applying  (1)  to  the  SBUV/2,  SAGE 
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II,  and  HALOE  ozone  profile  data  sets  (Figures  7  and  8).  Also  listed  is  the  mean  ozone 
concentration  (in  Dobson  Units,  DU,  per  km)  at  each  level  as  well  as  the  ENSO  ozone 
coefficients  calculated  by  Marsh  and  Garcia  for  the  1979-2003  period  (see  their  Figure 
4).  Although  direct  comparisons  are  not  straightforward  (due  to  the  different  regression 
models  that  were  employed),  at  several  pressure  levels  the  WACOM  coefficients  appear 
to  be  somewhat  larger  than  those  estimated  here  from  satellite  data  records.  Note  that 
the  largest  percentage  responses  are  found  in  the  lowermost  stratosphere  (~  70  hPa)  but 
the  ozone  concentration  at  this  level  is  relatively  low. 

Randel  et  al.  [2009b]  have  also  recently  estimated  SAGE  II  ENSO  ozone  regression 
coefficients  in  the  tropical  stratosphere  and  have  compared  these  coefficients  to  those 
calculated  from  WACOM  model  data.  They  employed  a  Multivariate  ENSO  Index  (MEI) 
rather  than  N3.4  and  the  SAGE  II  record  analyzed  was  slightly  longer  than  that  considered 
here.  (Note  that  the  mean  MEI  value  for  El-Nino  events  in  a  given  period  is  lower  than  the 
mean  N3.4  value  in  the  same  period.)  Just  above  the  equatorial  tropopause  (~  80  hPa; 
~  17.7  km),  a  large  percentage  response  maximum  was  obtained  with  a  peak  amplitude 
of  >  7  %/MEI  using  either  SAGE  data  or  WACOM  data  (compare  with  Figure  8a). 
However,  the  ozone  concentration  at  this  level  (~  2.3  DU/krn)  is  lower  than  at  higher 
altitudes  in  the  lower  stratosphere  (peaking  at  nearly  18  DU/krn  near  20  hPa;  see  Table  1). 
Similar  coherent  patterns  were  also  found  in  both  the  SAGE  II  and  WACOM  coefficients 
at  middle  latitudes  but  these  coherent  patterns  were  at  altitudes  below  ~  15  km  where 
the  mean  ozone  concentration  is  again  low.  Above  20  km  altitude  where  the  tropical 
ozone  concentration  is  large,  a  careful  reading  of  their  Figure  4  shows  that  the  ENSO 
ozone  coefficients  estimated  from  WACOM  model  data  are  somewhat  larger  than  those 
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calculated  from  SAGE  II  data.  For  example,  at  50  hPa,  the  equatorial  SAGE  II  ENSO 
coefficient  is  ~  2  %/MEI  while  that  estimated  from  WACCM  data  is  ~  3.4  %/MEI.  At 
40  hPa,  the  equatorial  SAGE  II  value  is  ~  1.1%/MEI  while  the  corresponding  WACCM 
value  is  nearly  2%/MEI.  At  30  hPa  where  the  tropical  ozone  concentration  approaches 
16  DU/km,  the  equatorial  SAGE  II  value  is  less  than  1%/MEI  while  the  corresponding 
WACCM  value  is  ~  2%/MEI.  Therefore,  at  levels  in  the  lower  stratosphere  where  the 
ozone  concentration  is  large,  their  results  appear  to  be  qualitatively  similar  to  those 
shown  in  Table  1. 

In  order  to  estimate  the  response  of  zonal  mean  stratospheric  temperature  and  zonal 
wind  to  ENSO,  we  consider  the  ECMWF  ERA-40  temperature  and  derived  zonal  wind 
data  set  extending  from  1979  through  2001.  Annual  mean  results  from  the  multiple  linear 
regression  analysis  are  shown  in  Figure  9.  Latitudes  from  60° S  to  60°N  are  given  to  be 
consistent  with  the  ozone  response  results  of  Figures  7  and  8.  Regression  coefficients 
are  shown  at  pressures  from  100  to  1  hPa,  which  differ  slightly  from  those  for  the  ozone 
response  plots  (50  to  0.5  hPa).  In  the  tropics,  the  ENSO  temperature  plot  (Figure  9a) 
is  characterized  by  a  negative  response  in  the  lowermost  stratosphere  (centered  in  the 
50  to  70  hPa  range).  This  negative  response  is  similar  to  the  lower  stratospheric  ENSO 
temperature  response  estimate  derived  originally  from  radiosonde  data  by  Reid  et  al. 
[1989]  and  Reid  [1994]  (see  also  Free  et  al.  [2009]  and  Randel  et  al.  [2009b]).  However, 
the  amplitude  of  the  radiosonde-derived  temperature  response  in  the  western  tropical 
Pacific  was  of  the  order  of  —1  K  for  a  change  of  one  standard  deviation  of  the  El-Nino 
SST  index  in  the  40  to  80  hPa  range  (comparable  to  a  change  in  N3.4  of  1  unit)  [Reid, 
1994],  In  contrast,  the  zonally  averaged  amplitude  derived  from  the  ERA-40  data  is  only 
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about  —0.1  K  for  an  increase  in  N3.4  of  1  unit.  The  zonal  wind  response  to  positive  ENSO 
forcing  (Figure  9b)  is  characterized  by  enhancements  of  the  subtropical  upper  tropospheric 
jets. 

The  annual  mean  TOMS/SBUV  total  ozone  ENSO  regression  coefficient  is  plotted  in 
Figure  10.  The  largest  region  of  statistically  significant  response  has  a  negative  amplitude 
and  is  located  in  the  tropical  eastern  Pacific  (~  30°S  to  30°N;  ~  80°W  to  180°W).  Note 
the  two  negative  anomaly  maxima  centered  on  the  equator  in  the  tropics.  This  pattern 
resembles  that  obtained  from  analyses  of  Microwave  Sounding  Unit  temperature  data 
[Yulaeva  and  Wallace,  1994;  Calvo-Fernandez  et  ah,  2004];  the  latter  authors  interpret 
this  structure  as  a  wave  pattern  that  mirrors  the  structure  in  the  troposphere.  Positive 
responses  with  amplitudes  comparable  to  those  of  the  tropical  negative  responses  are 
obtained  at  middle  latitudes  in  the  same  longitude  sector.  As  discussed  further  in  section 
5,  this  latitude  dependence  implies  ENSO-rclated  influences  on  the  residual  meridional 
circulation  (see  also  Free  et  al.  [2009],  and  section  5.1).  The  seasonal  dependence  of 
the  total  ozone  ENSO  coefficients  is  shown  in  Figure  11.  It  is  evident  from  these  plots 
that  the  negative  tropical  eastern  Pacific  response  is  most  pronounced  in  boreal  winter 
(DJF)  and  is  least  pronounced  in  summer  (JJA).  The  midlatitude  positive  responses  occur 
mainly  during  boreal  winter  and  spring.  In  particular,  during  boreal  winter,  the  strongest 
positive  maxima  in  the  Northern  Hemisphere  (up  to  2%  per  unit  of  N3.4)  occur  in  the 
northern  Pacific  region  south  of  Alaska.  These  northern  Pacific  response  patterns  may 
reflect  ENSO  influences  on  regional  stratospheric  circulation  such  as  the  strength  of  the 
Aleutian  stratospheric  anti-cyclone  in  addition  to  influences  on  the  residual  circulation. 
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4.  Model  Simulations  for  Warm/Cold  ENSO  Conditions 

Before  discussing  physical  mechanisms  in  section  5,  it  is  useful  to  compare  the  annual 
mean  ENSO  ozone  and  temperature  regression  coefficients  (Figures  7,  8,  and  9a)  with 
representative  model  simulations  for  strong/weak  ENSO  conditions.  For  this  purpose,  we 
use  the  Navy’s  Operational  Global  Atmospheric  Prediction  System  (NOGAPS),  which 
combines  a  global  spectral  forecast  model  with  a  3D  variational  data  assimilation  system. 
The  present  study  uses  a  high-altitude  version  of  NOGAPS,  NOGAPS-ALPHA  (Advanced 
Level  Physics-High  Altitude),  with  parameterized  chemistry  [McCormack  et  al.,  2006; 
Hoppcl  et  al.,  2008;  Eckermann  et  ah,  2008]. 

The  version  of  NOGAPS-ALPHA  used  here  has  68  (L68)  hybrid  a  —  p  vertical  levels 
extending  from  the  surface  to  5xl0-4  hPa  (~  100  km)  with  ~  2  km  spacing  throughout 
most  of  the  stratosphere  [Eckermann,  2009].  It  has  an  effective  Gaussian  horizontal  grid 
spacing  of  1.5°  in  latitude/longitude  (triangular  spectral  truncation  up  to  wave  number 
79).  Shortwave  heating  and  longwave  cooling  rates  due  to  O3,  CO2,  and  H20  from  the 
surface  to  the  lower  mesosphere  are  computed  using  the  method  of  Chou  et  al.  [2001] 
and  Chou  and  Suarez,  2002],  Above  90  km,  non-local  thermodynamic  equilibrium  (LTE) 
cooling  rates  due  to  C02  are  calculated  [Fomichev  et  al.,  1998].  Both  O3  and  H20  are 
prognostic  variables  in  the  forecast  model.  Production  and  loss  of  O3  and  H20  are  de¬ 
termined  using  a  linearized  photochemistry  parameterization  [McCormack  et  ah,  2006; 
2008].  The  model  is  therefore  not  a  fully  coupled  CCM.  In  particular,  while  photochemi¬ 
cal  effects  of  temperature  changes  are  approximately  accounted  for,  advection  of  long-lived 
trace  species  such  as  NOx,  which  significantly  influence  ozone  concentrations  in  the  middle 
and  upper  stratosphere,  is  not  accounted  for  in  the  model. 
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The  model’s  T79  spectral  truncation  cannot  explicitly  resolve  small-scale  gravity  wave 
breaking  in  the  middle  atmosphere,  which  exerts  a  drag  force  on  the  background  zonal 
winds  and  drives  the  circulation  away  from  a  thermal  equilibrium  state.  Therefore,  the 
model  uses  the  parameterizations  of  Webster  et  al.  [2003]  and  Garcia  et  al.  [2007]  to 
describe  the  effects  of  sub-grid  scale  orographic  and  non-orographic  gravity  wave  drag, 
respectively.  The  current  model  has  no  QBO,  either  forced  or  internally  generated. 

The  NOGAPS- ALPHA  forecast  model  uses  observed  or  climatological  distributions  of 
sea  surface  temperature  and  sea  ice  concentrations  as  lower  boundary  conditions.  To 
determine  the  response  of  NOGAPS- ALPHA  to  ENSO  forcing,  we  conducted  twelve  120- 
day  free-running  model  simulations  with  identical  initial  conditions.  The  first  ensemble 
of  six  simulations  used  operational  NOGAPS  SST  and  sea  ice  values  updated  every  12 
hours  for  the  period  2  December  2002  -  31  March  2003  (denoted  “DJF03”)  for  the  lower 
boundary  condition.  During  this  period,  SST’s  in  the  eastern  tropical  Pacific  were  much 
warmer  than  normal  (N3.4  ~  1-22),  characteristic  of  El-Nino  conditions  (see  Figure  12). 
The  second  ensemble  of  six  simulations  used  NOGAPS  SST  and  sea  ice  for  the  period 
2  December  2007  -  31  March  2008  (“DJF08”),  when  the  SST’s  were  cooler  than  normal 
(N3.4  ~  -1.72),  characteristic  of  La  Nina.  Initial  conditions  for  these  simulations  were 
generated  from  a  4-year  model  integration  using  monthly  mean  climatological  SST  and 
ice  fields  for  boundary  conditions  to  allow  the  model  to  reach  a  balanced  state.  Model 
output  was  saved  every  12  hours  throughout  the  simulations. 

Figures  13a, b,c  respectively  plot  the  zonal  mean  ozone  mixing  ratio  (ppmv),  tempera¬ 
ture  (Kelvin),  and  vertical  velocity  (nnn/s)  differences  (DJF03  —  DJF08)  as  a  function 
of  time  between  the  two  ensemble  means,  averaged  over  the  10°S  to  10°N  latitude  range. 
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Figure  14  replots  the  zonal  mean  ozone  mixing  ratio  differences  in  units  of  per  cent  to 

allow  more  direct  comparisons  to  the  observational  results  of  Figures  7  and  8  as  well  as  to 
the  SAGE  II  results  of  Randel  et  ah  [2009b],  Since  the  difference  in  N3.4  between  the  two 
time  periods  is  ~  3,  all  values  should  be  divided  by  3  to  estimate  the  model  response  to 
a  unit  change  in  N3.4.  In  the  last  column  of  Table  1,  we  list  the  ozone  differences  during 
the  final  60  days  of  the  integrations  averaged  over  24°S  to  24°N  and  expressed  in  units  of 
%/N3.4. 

In  the  case  of  the  ozone  difference  plots  (Figures  13a  and  14),  toward  the  end  of  the 
simulations,  ozone  mixing  ratio  differences  in  the  stratosphere  below  the  15  hPa  level 
are  generally  negative  and  range  from  ~  —0.13  pprnv  (~  2%)  near  20  hPa  to  ~  —0.06 
ppmv  (~  5%)  near  50  hPa,  i.e.,  the  mixing  ratios  are  reduced  under  El-Nino  conditions. 
These  differences  are  significant  at  >95%  confidence  as  verified  by  a  Student’s  t-test.  The 
corresponding  reductions  for  a  1-unit  decrease  in  N3.4  range  from  ~  0.04  to  0.02  ppmv.  As 
shown  in  Table  1,  the  mean  per  cent  differences  per  unit  of  N3.4  during  the  last  60  days  of 
the  model  integrations  averaged  over  24°S  to  24°N  range  from  ~  —0.5  %/N3.4  at  20  hPa 
to  ~  —1.2  %/N3.4  at  50  hPa.  These  model  sensitivities  are  comparable  in  magnitude  to 
the  observational  (e.g.,  SAGE  II)  ENSO  ozone  regression  coefficients  and  are  significantly 
less  than  the  WACCM  ENSO  coefficients  at  most  levels  in  the  lower  stratosphere  where 
the  mean  ozone  concentration  is  large. 

At  and  above  the  10  hPa  level,  mixing  ratio  differences  are  generally  positive;  mean 
differences  in  the  5-10  hPa  pressure  range  are  ~  +0.15  ppmv  (~  0.05  ppmv  for  a  1  unit 
change  in  N3.4),  implying  ozone  increases  under  El-Nino  conditions.  These  maximum  dif¬ 
ferences  are  again  significant  at  >  95%  conhdence.  The  model  response  vertical  structure 
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at  near-equatorial  latitudes  is  similar  to  that  of  the  observational  ENSO  ozone  regression 
coefficients  shown  in  Figures  7  and  8.  For  typical  climatological  mixing  ratios  at  5  -  10 
hPa  of  ~  9  ppmv,  the  ozone  increase  derived  from  SBUV/2  data  is  about  0.5%  or  ~  0.05 
ppmv,  not  much  less  than  the  mean  amplitudes  seen  in  Figures  13a  and  14  after  dividing 
by  3. 

In  the  case  of  the  temperature  difference  plot  (Figure  13b),  model  equatorial  temper¬ 
ature  differences  during  the  last  half  of  the  simulations  are  consistently  largest  in  the  50 
-  70  hPa  pressure  range  and  average  about  —2.5  K  (~  0.8  K  per  unit  of  N3.4).  These 
differences  are  significant  at  >  99%  confidence  according  to  a  Student’s  t-test.  As  seen 
in  Figure  9a,  the  observationally  derived  equatorial  temperature  difference  is  also  largest 
in  the  50  -  70  hPa  range  but  the  amplitude  is  only  —0.1  K  per  unit  of  N3.4.  The  model 
differences  also  show  lesser  decreases  of  as  much  as  1.5  K  (~  0.5  K  per  unit  of  N3.4) 
throughout  the  middle  and  upper  stratosphere,  which  are  not  seen  in  the  regression  re¬ 
sults.  On  the  other  hand,  as  also  noted  in  section  3,  Reid  [1994]  reported  temperature 
decreases  in  the  western  Pacific  lower  stratosphere  of  the  order  of  —  1  K  for  increases  of 
one  standard  deviation  of  an  SST  index  that  may  be  comparable  to  N3.4.  Also,  Randel 
et  al.  [2009b]  report  peak  equatorial  temperature  decreases  of  order  —0.5  K/MEI  in  the 
same  pressure  range  using  an  updated  and  adjusted  radiosonde  data  set. 

In  the  case  of  the  vertical  velocity  difference  plot  (Figure  13c),  toward  the  end  of  the 
simulations,  model  vertical  velocity  differences  are  generally  positive  throughout  most  of 
the  stratosphere  and  are  also  positive  in  the  upper  troposphere,  implying  net  tropical 
upwclling  during  warm  events  relative  to  cold  events.  Most  of  these  velocity  differences 
are  significant  at  >  95%  confidence  according  to  the  Student’s  t-test.  ft  is  unlikely  that 
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the  lack  of  significant  positive  velocities  in  the  lowermost  stratosphere  (~  60  to  100  hPa)  is 
related  in  any  way  to  the  model  formulation.  The  model  employs  a  single  hybrid  vertical 
coordinate  that  is  terrain  following  very  close  to  the  surface  and  isobaric  in  the  upper  levels 
[Eckermann,  2009].  Its  use  reduces  stratospheric  errors  in  divergence  fields  compared  to 
other  coordinate  systems,  and  so  is  well  suited  for  simulating  lower  stratospheric  vertical 
velocity  anomalies  related  to  higher  SST’s. 

It  should  be  emphasized  that  the  absence  of  a  QBO  in  the  model  undoubtedly  sim¬ 
plifies  the  model  response  compared  to  that  derived  from  observations.  In  addition,  the 
SST’s  drive  the  model  through  a  parameterization  of  cumulus  convection,  which  is  not 
straightforward  to  implement  in  such  a  model.  Nevertheless,  despite  these  caveats,  the 
approximate  agreement  of  the  NOGAPS- ALPHA  ENSO  ozone  sensitivities  with  observa¬ 
tions  (Table  1)  and  the  similar  altitude  dependences  of  the  modeled  and  observed  ENSO 
ozone  and  temperature  anomalies  (Figures  7,  8,  and  9a)  is  encouraging. 


5.  Physical  Interpretation  of  the  ENSO  Regression  Results 
5.1.  The  Lower  Stratospheric  Response 

As  seen  in  the  ERA-40  analysis  of  Figure  9a  and  in  the  NOGAPS-ALPHA  temperature 
differences  in  Figure  13b,  the  strongest  thermal  response  to  ENSO  forcing  in  the  tropical 
stratosphere  consists  of  a  pronounced  cooling  during  warm  events  centered  near  50  -  70 
hPa.  Reid  et  al.  [1989]  and  Reid  [1994]  first  reported  evidence  for  this  cooling  response 
using  radiosonde  data.  They  found  that  ENSO-related  cooling  effectively  amplified  the 
annual  temperature  cycle  in  this  region,  deepening  the  temperature  minimum  of  the  cycle, 
which  occurs  in  boreal  winter,  while  leaving  the  maximum  nearly  the  same.  Possible 
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mechanisms  discussed  at  that  time  for  producing  the  cooling  response  included:  (a)  the 
radiative  effect  of  optically  thick  cirrus  clouds  associated  with  enhanced  convective  activity 
during  warm  events  [Webster  and  Stephens,  1980];  and  (b)  variations  in  the  intensity  of 
the  tropospheric  Hadley  circulation  associated  with  increased  SST’s  [e.g.,  Bjerknes,  1966], 
which  could  produce  adiabatic  temperature  decreases  in  the  lowermost  stratosphere. 

As  shown  in  Figure  13c,  while  increased  upwclling  occurs  in  both  the  upper  troposphere 
and  in  the  lower  stratosphere  (50  hPa  and  above),  vertical  velocity  anomalies  tend  to  be 
quite  weak  near  60  hPa  where  the  cooling  maximizes  in  both  the  model  simulations 
(Figure  13b)  and  in  the  observations  (Figure  9a).  On  this  basis,  we  suggest  that  the 
strong  warm  event  cooling  seen  in  the  tropical  lowermost  stratosphere  is  probably  mostly 
radiative  in  origin,  deriving  from  decreases  in  upwelling  infrared  radiation  reaching  the 
lower  stratosphere  associated  with  increased  clouds  and  higher  cloud  tops  during  El-Nino 
events.  As  noted  by  Reid  [1994],  such  an  explanation  is  thermodynamically  reasonable 
since  the  tropospheric  heating  during  warm  events  must  be  compensated  by  a  cooling  at 
higher  altitudes  if  the  effects  of  ENSO  vanish  above  some  level. 

On  the  other  hand,  the  negative  ozone  ENSO  coefficients  in  the  lower  stratosphere  (Fig¬ 
ures  7  and  8)  and  the  mildly  reduced  temperatures  above  ~  50  hPa  seen  in  the  model 
simulations  of  Figure  13b  appear  to  be  mainly  caused  by  increased  upwelling  in  the  trop¬ 
ical  stratosphere  (Figure  13c).  Free  et  al.  [2009]  have  recently  found  evidence  that  the 
tropical  lower  stratospheric  cooling  in  boreal  winter  during  warm  events  is  accompanied 
by  a  strong  late  winter  warming  at  high  latitudes.  The  high-latitude  signal  has  charac¬ 
teristics  similar  to  those  of  sudden  stratospheric  warmings,  which  are  accompanied  by 
accelerations  of  the  Brewer-Dobson  circulation.  Consistently,  as  discussed  in  section  3, 
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column  ozone  ENSO  regression  coefficients  are  negative  at  low  latitudes  and  positive  at 
high  latitudes  in  the  eastern  Pacific  region  (Figure  10).  Since  the  Brewer-Dobson  circula¬ 
tion  is  driven  primarily  by  planetary  wave  activity,  it  can  be  inferred  that  the  enhanced 
tropical  upwelling  is  a  result  of  enhanced  Rossby  wave  activity  associated  with  the  ENSO- 
induced  tropospheric  circulation  changes.  In  support  of  this  explanation,  several  recent 
modeling  studies  have  shown  that  increased  SST’s  can  modify  quasi-stationary  planetary 
waves  in  the  tropics  with  effects  on  the  residual  circulation  [e.g.,  Deckert  and  Dameris, 
2008].  In  particular,  Garcia-Herrera  et  al.  [2006]  have  studied  the  propagation  of  tro¬ 
pospheric  ENSO-related  planetary  wave  disturbances  into  the  winter  stratosphere  using 
both  general  circulation  models  and  the  ERA-40  data.  In  general,  warm  ENSO  events 
tend  to  maximize  in  boreal  winter  and  produce  zonally  asymmetric  tropospheric  tem¬ 
perature  anomalies  in  the  tropics  that  propagate  in  the  troposphere  to  the  extratropics. 
The  resulting  midlatitude  temperature  anomalies  then  propagate  vertically  as  planetary 
Rossby  waves  into  the  midlatitude  stratosphere  where  they  decelerate  the  westerly  flow, 
inducing  an  acceleration  of  the  residual  circulation.  The  net  consequence  for  the  tropical 
stratosphere  is  increased  upwelling  rates.  This  mechanism  is  currently  being  examined 
further  using  the  NOGAPS-ALPHA  simulations  described  here. 

5.2.  The  Middle  Stratospheric  Response 

It  remains  to  consider  the  origin  of  the  positive  tropical  ENSO  ozone  response  in  the 
middle  stratosphere  (5  to  10  hPa)  seen  in  the  regression  results  of  Figures  7  and  8  and 
in  the  model  results  of  Figure  13a.  A  similar  positive  ENSO  ozone  response  centered 
near  10  hPa  was  obtained  in  the  analysis  of  WACOM  model  data  by  Marsh  and  Garcia 
[2007]  who  attributed  it  to  the  photochemical  effects  of  adiabatic  temperature  changes 
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associated  with  increased  tropical  upwelling  after  warm  events.  As  was  emphasized  in 
the  model  description  of  section  4,  the  NOGAPS-ALPHA  model  currently  has  only  a 
parameterized  chemistry  scheme  and  therefore  accounts  only  for  ozone  changes  resulting 
from  dynamical  transport  effects  on  temperature,  which  in  turn  affect  reaction  rates  that 
determine  the  ozone  photochemical  balance.  The  positive  ozone  changes  seen  in  Figure  13a 
above  10  hPa  are  therefore  caused  primarily  by  adiabatic  temperature  decreases  (Figure 
13b)  associated  with  the  upwelling  rate  increases  during  warm  events  (Figure  13c). 

However,  in  the  stratosphere  (and  in  full  CCM’s  such  as  WACOM),  changes  in  upwelling 
rates  near  10  hPa  produce  transport-induced  changes  in  other  trace  gases  with  longer 
lifetimes  than  ozone.  In  particular,  transport-induced  changes  in  odd  nitrogen  occur, 
which  can  strongly  influence  the  ozone  photochemical  balance  by  modifying  catalytic  loss 
rates.  To  illustrate  this,  Figure  15  compares  zonal  averages  of  NOx  and  O3  between  5°S 
and  5°N  at  10  hPa  as  calculated  from  UARS  HALOE  sunset  data  over  a  12-year  period. 
A  clear  NOx  QBO  is  evident  that  is  approximately  out  of  phase  with  (and  significantly 
correlated  with)  the  ozone  QBO  (R  =  -0.94)  [Randel  and  Wu,  1996;  Baldwin  et  ah,  2001]. 
It  is  therefore  likely  that  increased  tropical  upwelling  following  warm  ENSO  events  would 
lead  to  reduced  NOx  concentrations  and  increased  ozone  in  the  middle  stratosphere.  This 
mechanism  may  even  dominate  over  the  photochemical  effects  of  adiabatic  temperature 
changes,  as  is  apparently  true  for  the  QBO  [Chipperfield  et  ah,  1994], 

5.3.  Implications  for  the  Solar  Cycle  Response 

It  has  previously  been  inferred  from  observations  that  increased  tropical  upwelling  oc¬ 
curs  near  or  approaching  solar  minima  [Kodera  and  Kuroda,  2002;  Hood  and  Soukharev, 
2003].  In  the  tropical  lower  stratosphere,  such  increased  upwelling  would  reduce  ozone 
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mixing  ratios  near  solar  minima,  consistent  with  the  observed  second  lower  stratospheric 
ozone  response  maximum  (e.g.,  Figure  3).  In  the  tropical  middle  stratosphere,  based  on 
the  interpretations  given  in  the  previous  subsection,  it  is  possible  that  the  absence  of  a 
detectable  solar  cycle  ozone  response  is  also  due  to  increased  upwclling  near  solar  minima, 
which  would  advectively  decrease  both  the  temperature  and  the  odd  nitrogen  abundance, 
leading  to  photochemical  ozone  increases.  Such  an  ozone  increase  near  solar  minima  would 
tend  to  cancel  the  ozone  decrease  at  that  time  caused  by  reduced  photolysis  of  molecular 
oxygen  by  solar  UV  radiation.  The  net  result  would  be  an  insignificant  ozone  response  in 
the  tropical  middle  stratosphere. 

To  examine  this  hypothesis  further,  it  is  useful  to  consider  again  the  odd  nitrogen  and 
ozone  time  series  shown  in  Figure  15  at  the  location  of  the  ozone  response  minimum  (~  10 
hPa,  equator).  According  to  these  data,  the  interannual  variability  of  both  equatorial  odd 
nitrogen  and  ozone  at  this  location  is  dominated  by  an  irregular  QBO,  making  detection 
of  the  solar  cycle  component  of  variability  more  difficult.  Model  calculations  (see,  e.g., 
Figure  14  of  SH06)  indicate  that  the  photochemical  ozone  increase  from  solar  minimum  to 
maximum  at  this  level  is  only  ~  2%.  In  contrast,  the  peak-to-peak  ozone  QBO  amplitude 
seen  in  Figure  15b  is  in  the  range  of  10  to  15%.  Regressing  the  NOx  time  series  of  Figure 
15a  against  the  ozone  time  series  of  Figure  15b  yields  ANOx  ~  — 2.9AO3.  Consequently, 
a  transport-induced  increase  in  NOx  of  ~  5.8%  from  solar  minimum  to  maximum  may  be 
sufficient  to  cancel  a  direct  solar  UV-induced  ozone  increase  of  2%  at  that  location.  Such 
an  NOx  increase  is  much  less  than  the  peak-to-peak  QBO-induced  amplitude  (>  40%) 
seen  in  Figure  15a  and  is  therefore  difficult  to  rule  out.  For  example,  simple  averaging 
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of  the  data  of  Figure  15a  for  various  time  periods  representative  of  solar  minimum  and 
maximum  conditions  yields  N0X  increases  in  the  range  of  1  to  5%. 

6.  Summary  and  Discussion 

The  analysis  of  NCAR  WACOM  model  data  for  the  1979-2003  time  period  by  Marsh 
and  Garcia  [2007]  was  valuable  for  drawing  attention  to  possible  aliasing  of  solar  ozone  re¬ 
gression  coefficients  by  ENSO.  In  particular,  they  found  that  inclusion  of  an  ENSO  term 
in  a  regression  model  can  produce  large  (~  35%)  reductions  in  WACOM  annual  mean 
lower  stratospheric  ozone  solar  regression  coefficients  (see  their  Figure  4).  However,  the 
present  results  indicate  that  when  observed  ozone  records  (SBUV/2,  SAGE  II,  HALOE) 
are  analyzed  using  a  similar  statistical  model,  only  minor  changes  to  annual  mean  coef¬ 
ficients  (Figures  3  and  4)  are  obtained  when  an  ENSO  term  is  included.  The  same  is 
true  when  seasonal  data  are  analyzed  for  three  of  the  four  seasons.  Only  when  data  are 
analyzed  for  boreal  winter  (DJF),  when  ENSO  events  tend  to  peak,  is  evidence  obtained 
for  limited  (<  25%)  ENSO  aliasing  of  solar  regression  coefficients  within  several  latitude, 
longitude,  and  pressure  ranges  (Figures  5  and  6). 

As  noted  in  section  2.2,  a  partial  explanation  for  the  stronger  influence  of  ENSO  on  solar 
regression  coefficients  derived  from  WACOM  model  data  is  the  presence  in  observational 
records  of  other  sources  of  interannual  variability,  e.g.,  an  irregular  QBO,  which  would 
tend  to  obscure  the  ENSO  signal.  In  addition,  as  discussed  in  section  3,  it  is  likely  that  the 
ENSO  signal  in  WACOM  is  somewhat  stronger  than  is  seen  in  observations  at  a  number 
of  levels  in  the  lower  stratosphere  where  the  ozone  concentration  is  high  (see  Table  1  and 
Figure  4  of  Randel  et  al.  [2009b]).  Although  the  latter  authors  have  shown  that  similar 
coherent  positive  and  negative  patterns  of  per  cent  ozone  ENSO  response  are  seen  in  both 
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WACOM  and  in  SAGE  II  data,  most  of  the  coherent  patterns  occur  at  altitudes  of  18 
km  and  below  where  the  ozone  concentration  is  low.  Consistent  with  these  inferences, 
WACOM  model  ozone  time  series  near  50  hPa  (Figure  2  of  Marsh  and  Garcia)  correlate 
strongly  with  N3.4  (R  =  —0.7  at  4  months  lag)  while  the  ENSO  signal  in  observed  ozone 
records  in  the  tropics  is  weaker  (e.g.,  R  =  —0.37  at  4  months  lag  for  the  column  ozone 
time  series  in  Figure  2). 

These  results  do  not  imply  that  the  estimated  solar  and  ENSO  regression  coefficients 
are  highly  accurate.  The  accuracy  of  the  solar  coefficients,  in  particular,  depends  also  on 
the  length  of  the  time  series,  which  for  current  satellite  datasets,  is  limited  to  no  more 
than  2  or  3  solar  cycles.  A  completely  accurate  determination  of  the  linear  response  may 
require  as  many  as  10  cycles  of  high-quality  global  measurements.  Further,  it  is  well 
known  that  the  QBO,  ENSO,  and  solar  cycle  responses  interact  and  their  effects  are  not 
linear,  especially  at  high  latitudes  [e.g.,  Garhnkel  and  Hartmann,  2007;  Calvo  et  al.,  2009; 
Labitze,  2004],  which  is  inconsistent  with  the  assumptions  of  the  multiple  linear  regression 
method.  Hence,  even  with  an  infinitely  long  time  series,  the  solar  and  ENSO  regression 
coefficients  may  not  attain  complete  accuracy. 

As  found  in  section  3,  the  ENSO  ozone  regression  coefficient  derived  from  SBUV/2, 
SAGE  II,  HALOE,  and  TOMS  data  is  characterized  by  a  negative  response  in  the  tropical 
lower  stratosphere,  mainly  in  the  eastern  Pacific  region  over  a  wide  range  of  altitudes 
(Figures  7,  8,  and  10).  The  temperature  coefficient  is  characterized  by  a  negative  response 
in  the  lower  stratosphere  that  is  statistically  significant  within  a  narrow  pressure  range 
centered  at  50  —  70  hPa.  The  zonal  wind  response  consists  mainly  of  an  acceleration  of 
the  equatorward  flanks  of  the  tropospheric  subtropical  jets  (Figure  9).  The  difference 
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between  NOGAPS-ALPHA  model  runs  for  continuous  warm  and  cold  ENSO  conditions 
at  low  latitudes  also  has  characteristics  that  are  very  similar  to  the  calculated  tropical 
ENSO  ozone  and  temperature  regression  coefficients  (section  4;  Figure  13;  Table  1). 

As  discussed  in  section  5.1,  when  both  the  observational  and  model  results  are  con¬ 
sidered  together,  the  pronounced  cooling  at  50  —  70  hPa  during  warm  events  can  be 
interpreted  as  primarily  a  radiative  response  to  increased  cloudiness  in  the  upper  tropo¬ 
sphere.  However,  the  negative  ozone  response  and  weaker  cooling  observed  throughout 
the  tropical  lower  stratosphere  are  attributable  to  an  indirect  dynamical  effect,  i.e.,  an 
acceleration  of  the  Brewer-Dobson  circulation,  which  is  a  result  of  enhanced  Rossby  wave 
activity  associated  with  the  tropical  tropospheric  circulation  changes  (see,  e.g.,  Yulaeva  et 
ah,  1994;  Garcia-Herrera  et  ah,  2006;  Deckert  and  Dameris,  2008;  Randel  et  ah,  2008].  As 
noted  previously  by  Shibata  and  Deushi  [2008],  the  observed  zonal  wind  ENSO  response 
can  be  interpreted  as  essentially  a  thermal  wind  response  to  increased  latent  heat  release 
in  the  tropical  troposphere.  As  discussed  in  section  5.2,  the  positive  ozone  response  ob¬ 
served  in  the  tropical  middle  stratosphere  (centered  near  10  hPa)  is,  by  analogy  to  the 
ozone  QBO,  attributable  to  advective  temperature  and  odd  nitrogen  decreases,  which 
would  lead  to  photochemical  ozone  increases. 

In  section  5.3,  it  was  first  noted  that  previous  work  has  indicated  that  the  tropical 
upwclling  rate  increases  near  or  approaching  solar  minima  [e.g.,  Kodera  and  Kuroda, 
2002;  Hood  and  Soukharev,  2003].  Such  upwclling  rate  changes  would  lead  to  transport- 
induced  ozone  increases  in  the  lower  stratosphere  near  solar  maxima,  consistent  with  the 
observed  total  ozone  variation  (e.g.,  Figure  lb).  It  was  then  suggested  that  the  observed 
insignificant  ozone  response  to  the  solar  cycle  in  the  tropical  middle  stratosphere  could 
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be  a  consequence  of  transport-induced  increases  in  temperature  and  odd  nitrogen  mixing 
ratios  (by  analogy  to  those  that  occur  in  cold  ENSO  events).  The  latter  would  catalytically 
decrease  the  ozone  mixing  ratio,  tending  to  cancel  out  the  ozone  increase  resulting  from 
direct  solar  UV-induced  increases  in  ozone  production. 
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Figure  1.  Comparisons  of  the  solar  Mg  II  UV  index  with  area- weighted  averages  over  two 
latitude  ranges  of  the  monthly  mean  version  8  TOMS/SBUV(/2)  total  ozone  anomaly.  Also 
indicated  are  the  times  of  the  El  Chichon  (C)  and  Pinatubo  (P)  volcanic  eruptions. 

Figure  2.  Comparison  of  the  Nino  3.4  index  with  the  TOMS/SBUV(/2)  total  ozone  anomaly 
averaged  over  tropical  latitudes.  Also  indicated  are  the  correlation  coefficients  (R)  at  zero  lag 
and  at  4  months  lag  (the  lag  at  which  \R\  maximizes). 

Figure  3.  Comparison  of  annual  mean  solar  cycle  ozone  regression  coefficients  calculated  from 
SBUV(/2)  data  (a)  using  the  regression  model  (1)  without  an  ENSO  term  and  (b)  using  the 
regression  model  (1)  with  an  ENSO  term  at  optimum  lag.  The  contour  interval  is  0.5%  and 
shaded  areas  are  significant  at  95%  conhdence. 

Figure  4.  Comparison  of  tropically  averaged  solar  cycle  ozone  regression  coefficients  calculated 
from  the  three  satellite  data  sets  analyzed  by  Soukharev  and  Hood  [2006]  (top  panel)  using  the 
regression  model  (1)  without  an  ENSO  term  and  (bottom  panel)  using  the  regression  model  (1) 
with  an  ENSO  term  at  optimum  lag. 

Figure  5.  Same  format  as  Figure  3  but  for  boreal  winter  (December,  January,  February)  solar 
cycle  ozone  regression  coefficients. 

Figure  6.  Comparison  of  solar  cycle  column  ozone  regression  coefficients  for  boreal  winter 
(DJF)  calculated  from  Version  8  TOMS/SBUV  data  (a)  using  the  regression  model  (1)  without 
an  ENSO  term  and  (b)  using  the  model  (1)  with  an  ENSO  term  at  optimum  lag. 

Figure  7.  Per  cent  change  in  version  8  SBUV(/2)  ozone  for  a  1  unit  change  of  the  Nino  3.4 
index.  Shaded  areas  are  significant  at  the  95%  conhdence  level. 
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Figure  8.  Same  format  as  Figure  7  but  for  (a)  SAGE  ff  ozone  profile  data  over  the  1985-2003 
period  and  (b)  UARS  ffALOE  ozone  profile  data  over  the  1992-2003  period. 

Figure  9.  Annual  mean  ENSO  regression  coefficients  from  ERA-40  data  for  (a)  temperature 
in  Kelvin  per  unit  of  N3.4;  and  (b)  zonal  mean  zonal  wind  in  m/s  per  unit  of  N3.4.  Shaded  areas 
are  significant  at  the  95%  confidence  level. 

Figure  10.  Same  format  as  Figure  7  but  for  annual  mean  version  8  TOMS/SBUV(/2)  column 
ozone  data  as  a  function  of  latitude  and  longitude. 

Figure  11.  Seasonal  dependence  of  the  TOMS/SBUV(/2)  total  ozone  ENSO  regression 
coefficient.  From  top  to  bottom:  Winter  (DJF),  Spring  (MAM),  Summer  (JJA),  Fall  (SON). 
Shaded  areas  are  significant  at  95%  confidence  level. 

Figure  12.  Comparison  of  the  Nino  3.4  ENSO  index  for  the  two  northern  winter  periods  for 
which  simulations  were  performed  using  the  NOGAPS-ALPHA  model. 

Figure  13.  Difference  between  the  NOGAPS-ALPHA  ensemble  model  runs  for  warm  (DJF03) 
and  cold  (DJF08)  ENSO  conditions  as  a  function  of  time  (days),  averaged  over  the  10°S  to 
10°N  latitude  range,  for  (a)  ozone  mixing  ratio  (C.l.  =  0.05  ppmv);  (b)  temperature  (C.I.  =  0.5 
Kelvin);  and  (c)  zonally  averaged  vertical  velocity  (contours  at  ±0.1,  ±0.2,  ±0.4,  ±0.8,  ±1.2  ... 
mrn/s).  Blue  and  red  colors  indicate  negative  and  positive  values,  respectively.  All  values  should 
be  divided  by  ~  3  to  convert  to  units  per  unit  of  N3.4. 

Figure  14.  Same  as  Figure  13a  but  plotted  in  units  of  per  cent  (contours  at  ±1,  ±2,  ±3,  ±4, 
±5,  ±10,  and  ±20%).  Values  should  again  be  divided  by  ~  3  to  convert  to  per  cent  per  unit  of 
N3.4. 


DRAFT 


February  4,  2010,  10:02am 


DRAFT 


HOOD  ET  AL.:  DECADAL  VARIABILITY  OF  THE  TROPICAL  STRATOSPHERE 


X  -  45 


Figure  15.  Comparison  of  3-month  zonal  averages  of  UARS  HALOE  sunset  measurements  at 
10  hPa,  5°S  to  5°N,  of  (a)  NO  +  N02;  and  (b)  03  over  the  1992  to  2003  period. 
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Table  1.  Mean  O3  Concentration  and  ENSO  O3  Regression 
Coefficients,  24°S  -  24°N,  1979-2003 


p  [03]* 
hPa  DU /km 

SBUV 

%/N3.4 

SAGE  II 
%/N3.4 

HALOE 

%/N3.4 

WACCM** 

%/N3.4 

NOGAPST 

%/N3.4 

20 

17.7 

-0.69 

-0.37 

-0.68 

-0.5 

-0.5 

30 

15.8 

-0.86 

-0.50 

-0.72 

-1.5 

-0.6 

40 

13.0 

-0.85 

-0.57 

-0.47 

-2.2 

-0.8 

50 

9.6 

-0.75 

-1.25 

-0.16 

-2.8 

-1.2 

70 

3.4 

— 

-3.32 

-1.63 

-3.3 

-4.2 

*  Approximate  long-term  means  from  SAGE  II  data 
**  From  Figure  4  of  Marsh  and  Garcia  [2007] 

f  Estimated  from  mean  differences  over  the  last  60  days  of  the  integrations 
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